#include <limits>

double sum(double a, double b, double & roundoff);
void split(double a, std::size_t s, double & hi, double & lo);
double mul(double a, double b, double & roundoff);

